Exponential sensitivity to dephasing of electrical conduction through a quantum dot 
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According to random-matrix theory, interference effects in the conductance of a ballistic chaotic 
quantum dot should vanish oc (t$/td) p when the dephasing time becomes small compared to 
the mean dwell time td- Aleiner and Larkin have predicted that the power law crosses over to an 
exponential suppression oc exp(— te/t$) when drops below the Ehrenfest time te- We report the 
first observation of this crossover in a computer simulation of universal conductance fluctuations. 
Their theory also predicts an exponential suppression oc exp(— te/td) in the absence of dephasing 
— which is not observed. We show that the effective random-matrix theory proposed previously for 
quantum dots without dephasing explains both observations. 
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An instructive way to classify quantum interference ef- 
fects in mesoscopic conductors is to ask whether they 
depend exponentially or algebraically on the dephasing 
time T^. The Aharonov-Bohm effect is of the former 
class, while weak localization (WL) and universal con- 
ductance fluctuations (UCF) are of the latter class 0, • 
It is easy enough to understand the difference: On the 
one hand, Aharonov-Bohm oscillations in the magneto- 
conductance of a ring require phase coherence for a cer- 
tain minimal time i m i n (the time it takes to circulate once 
along the ring), which becomes exponentially improbable 
if < t m i n . On the other hand, WL and UCF in a dis- 
ordered quantum dot originate from multiple scattering 
on a broad range of time scales, not limited from below, 
and the superposition of exponents with a range of decay 
rates amounts to a power law decay. 

In a seminal paper |3|, Aleiner and Larkin have ar- 
gued that ballistic chaotic quantum dots are in a class of 
their own. In these systems the Ehrenfest time te intro- 
duces a lower limiting time scale for the interference ef- 
fects, which are exponentially suppressed if < te- The 
physical picture is that electron wave packets in a chaotic 
system can be described by a single classical trajectory 
for a time up to te 0- Both WL and UCF, however, re- 
quire that a wave packet splits into partial waves which 
follow different trajectories before interfering. Only the 
fraction exp(— te/t$) of electrons which have not yet de- 
phased at time te can therefore contribute to WL and 
UCF. 

The WL correction AG = (G) — G c \ is the deviation 
of the ensemble averaged conductance (G) (in zero mag- 
netic field) from the classical value G c \ = N/2. (We 
measure conductances in units of 2e 2 /h and assume an 
equal number of modes N 3> 1 in the two leads that con- 
nect the quantum dot to electron reservoirs.) The WL 
correction according to random-matrix theory (RMT), 
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has a power 



AGrmt = -j(1+td/t*) 
law suppression oc t^/td when becomes 



smaller than the mean dwell time td in the quantum dot 
0. Similarly, RMT predicts for the UCF a power law 
suppression cx (t^/td) 2 of the mean-squared sample-to- 
sample conductance fluctuations 0, , 
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with j3 = 2 (1) in the presence (absence) of a time- 
reversal-symmetry-breaking magnetic field. 

Aleiner and Larkin have calculated the r^-dependence 
of the WL correction, with the result 



AG = e- TE/T *e~ 2TE/TD AG 



RMT- 



(3) 



The two exponential suppression factors in Eq. J3J result 
from the absence of interfering trajectories for times be- 
low te- The first factor exp(— te/t^) accounts for the 
loss by dephasing and the second factor exp(— 2te/te) 
accounts for the loss by escape into one of the two leads. 
The UCF are expected to be suppressed similarly. 

The physical picture presented by Aleiner and Larkin 
is simple and supported by two independent calculations 
00. And yet, it has been questioned as a result of some 
very recent computer simulations of UCF 0, and WL 
in the absence of dephasing. The expected expo- 
nential reduction of quantum interference effects due to 
escape into the leads was not observed. In fact, both WL 
and UCF were found to be completely independent of te, 
even though the simulations extended to system sizes for 
which te was well above td- To explain these negative 
results, Jacquod and Sukhorukov invoked the effective 
RMT of Silvestrov et al. [llj . In that approach the elec- 
trons with dwell times > te are described by RMT with 
an effective number N c s — iVexp(— te/td) of modes. 
Then no r^-dependence is expected as long as iV c ff 3> 1 
— even if te 3> td- 

Since the predicted exponential reduction factor due 
to escape into the leads has not appeared in the simula- 
tions, it is natural to ask about the factor oc exp(— te/t^,) 
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due to dephasing. Does it exist? An experimental study 
of two-dimensional (2D) weak localization has concluded 
that it does |l2j, but since leads play no role in 2D 
these experiments can not really resolve the issue. In 
the absence of experiments on the zero-dimensional ge- 
ometry of a quantum dot, we have used computer sim- 
ulations to provide an answer. We find that a relatively 
small amount of dephasing is sufficient to introduce a 
marked r^-dependence of the UCF. Our observation can 
be explained by incorporating dephasing into the effec- 
tive RMT. We find that 



AG = e~ TE / T * AGrmt, 
Var G = e - 2TE/T * Var G RM t, 



(4) 
(5) 



and show that Eq. (jSJ provides a fitting-parameter-free 
description of the numerical data. 

We have introduced a dephasing lead 0] in the kicked 
rotator, which is the same dynamical system studied 
in Refs. H, 0, 0, Q in the absence of dephasing. 
The kicked rotator provides a stroboscopic description 
of chaotic scattering in a quantum dot |l5j , in the sense 
that the wave function is determined only at times which 
are multiples of a time To (which we set to unity). The 
mean dwell time td = M /2N = tt/N5 is the ratio of the 
dimension M of the Floquet matrix (corresponding to a 
mean level spacing S — 2n/M) and the dimension 27V of 
the scattering matrix (without the dephasing lead). The 
kicking strength K — 7.5 determines the Lyapunov ex- 
ponent A = \n(K/2) — 1.32. The Ehrenfest time is given 

by BG3 



TE 



X- 1 ln(jV 2 /M) if N > 







if N < 
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The dephasing lead increases the dimension of the scat- 
tering matrix S to M x M. It has the block form 
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FIG. 1: Variance of the conductance fluctuations, normalized 
by the RMT value as a function of the dimension M of the 
scattering matrix of the kicked rotator with a dephasing lead. 
Each data set is for a fixed value of the dwell time td = M /2N 
and dephasing time = T _1 (l — 2JV/M) -1 . The Lyapunov 
exponent A = 1.32 is kept the same for all data sets. The 
curves show the Ehrenfest time dependence iJKJ predicted by 
the effective RMT, without any fit parameter. 
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where the subscripts 1,2 label the two real leads and 
labels the dephasing lead. The two real iV-mode leads are 
coupled ballistically to the system, while the remaining 
M — 2N modes are coupled via a tunnel barrier. The 
dephasing rate 1/r^ = T(l — 2N/M) is proportional to 
the tunnel probability Y per mode. The dephasing lead 
is connected to an electron reservoir at a voltage which is 
adjusted so that no current is drawn. The conductance 
G is then determined by the coefficients Gy — ' • • - 
through Buttiker's formula [l3|. 



G = Gi 
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2N — Gn ~ G12 — G22 — G21 
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For r <C I the dephasing lead model is equivalent to 
the imaginary energy model of dephasing 5] , which is the 
model used by Aleiner and Larkin Q . (We will also make 
use of this equivalent representation later on.) There 
exist other models of dephasing in quantum transport 
[lsL ITgj . but for a comparison with Ref. Q our choice 
seems most appropriate. 

Since we need a relatively small Lyapunov exponent in 
order to reach a large enough Ehrenfest time, our sim- 
ulations are sensitive to short non-ergodic trajectories. 
These introduce an undesired dependence of the data on 
the position of the leads. Preliminary investigations in- 
dicated that UCF in a magnetic field was least sensitive 
to the lead positions, so we concentrate on that quantum 
interference effect in the numerics. The variance Var G of 
the conductance was calculated in an ensemble created 
by sampling 40 values of the quasi-energy. To determine 
the dependence on the lead positions we repeated the cal- 
culation for 40 different configurations of the leads. Error 
bars in the plots give the spread of the data. 
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FIG. 2: Pictorial representation of the effective RMT of a 
ballistic chaotic quantum dot. The part of phase space with 
dwell times > te is represented by a fictitious chaotic cavity 
(mean level spacing <5 c ff), connected to electron reservoirs by 
two long leads (N c ff propagating modes, one-way delay time 
te /2 for each mode) . The effective parameters are determined 
by N e fi/N — S/8 c g — exp(— te/td)- The scattering matrix 
of lead plus cavity is exp(iETE/h)SnMT(E), with 5'rmt(-B) 
distributed according to RMT. A finite dephasing time is 
introduced by the substitution E — > E + i%/2r^. The part of 
phase space with dwell times < te has a classical scattering 
matrix, which does not contribute to quantum interference 
effects. 



There are four time scales in the problem: A -1 , td, t$, 
and te- To isolate the te dependence we increase both M 
and N at constant ratio M/N and fixed K, T. Then only 
te varies. Results are shown in Fig. The variance of 
the conductance is divided by the RMT prediction 
with (3 = 2 because of broken time-reversal symmetry 
|2fij . We see that for 3> te there is no systematic de- 
pendence of UCF on the Ehrenfest time, consistent with 
Refs. However, an unambiguous ^-dependence 

appears for < te, regardless of whether is smaller 
or larger than tb- 

To explain the data in Fig. ^ we introduce dephas- 
ing into the effective RMT. For that purpose it is more 
convenient to use an imaginary energy than a dephasing 
lead, so we first make the connection between these two 
equivalent representations. There exists an exact corre- 
spondence for any TV which requires a re-injection 
step to ensure current conservation. For the case N ^> 1 
of interest here there is a simpler way. 

The coefficients G y = Gfj + G?. in Eq. © consist 
of a classical contribution G^J of order N plus a (sam- 
ple specific) quantum correction G?- of order unity. The 
classical contribution is 

Gf j = ±N(l+T D /T 4> )-\ for *,je {1,2}. (9) 

Substitution into Eq. JSJ) gives a classical conductance 
G c i = N/2 independent of dephasing — as it should be. 
To leading order in N we obtain the quantum correction 
to the conductance, 



= 4(^12 + G21 — Gn — G2 



(10) 



(Notice that the classical contribution drops out of the 
right-hand-side.) For f < 1 the effect of the dephasing 



lead on the coefficients Gy is equivalent to the addition 
of an imaginary part ih/2r^ to the energy. With the 
help of Eq. (fTTTf) we can compute the effect of dephasing 
on WL and UCF, 

AG = (G q (E + ih/2T c f,)), (11) 
VarG= ([G q (E + th/2T^)} 2 ) - (AG) 2 , (12) 

by averaging the scattering matrix at a complex energy 
without having to enforce current conservation. 

Effective RMT [Tl| is a phenomenological decompo- 
sition of the scattering matrix S(t) in the time domain 
into a classical deterministic part S c \ for t < te and a 
quantum part S q with RMT statistics for t > te, 



cm _ / s c\(t) if t <T E , 

" ,[T > \ S q (t) = Sn MT (t - t e ) if t>T E . 



(13) 



The RMT part S q couples to a reduced number N e g = 
N exp(— te I Try) of channels in each lead. The mean dwell 
time in the quantum dot of these channels is te+td- The 
classical part S c i couples to the remaining 2 (TV — TV c ff) 
channels. (See Ref. [TtI for an explicit construction of 

fid.) 

Only S q contributes to G q . Fourier transformation to 
the energy domain gives 



S q (E) 



TE/h S R MT(E), 



(14) 



where we have used that SrmtW = if t < 0. The ma- 
trix Srmt has the RMT statistics of a fictitious chaotic 
cavity with zero Ehrenfest time, 7V e g modes in each 
lead, and the same mean dwell time te as the real 
cavity (see Fig. For real energy the phase factor 
exp(iETE/K) is irrelevant, hence all r^-dependence is 
hidden in N e s and 5 e Q. Since AG and Var G are indepen- 
dent of these two parameters, they are also independent 
of te- The imaginary part ih/2T ( / > of the energy that 
represents the dephasing introduces a r^-dependence of 
G q oc exp(— te/t$). Insertion of this factor into Eqs. 
(jllfl and (|12|) yields the results |@} and JSJ given in the 
introduction. 

The curves in Fig.Ufollow from Eq. JSJl. They describe 
the simulation quite well — without any fit parameter. 
To test the agreement between simulation and effective 
RMT in a different way, we have collected all our data 
in Fig.|3]in a plot of — (r^/2) ln(VarG/VarGRMT) versus 
ln(7V 2 /M). According to Eq. JSJ this should be a plot 
of te versus ln(7V 2 /M), which in view of Eq. 10 is a 
straight line with slope 1/A = 0.76. There is considerable 
scatter of the data in Fig. [21 but the systematic parameter 
dependence is consistently described by the theory as N 
and M vary over two orders of magnitude. 

In conclusion, our findings explain the puzzling differ- 
ence in the outcome of previous experimental [l2j and 
numerical 0, 0, searches for the Ehrenfest time de- 
pendence of quantum interference effects in chaotic sys- 
tems: The experiments found a dependence while the 
computer simulations found none. We have identified the 
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FIG. 3: Four data sets of fixed t$, td, each consisting of 
a range of M between 10 2 and 2 ■ 10 4 , plotted on a double- 
logarithmic scale. The solid line with slope 1/A = 0.76 is the 
scaling predicted by Eqs. @ and |0. 

absence of dephasing in the simulations as the origin of 
the difference. By introducing dephasing into the simula- 
tion we recover the exponential te / ">> suppression factor 



predicted by Aleiner and Larkin 0. The effective RMT 
explains why this suppression factor is observed while the 
exponential te/td suppression factor of Eq. J3J) is not. 

It remains an outstanding theoretical challenge to pro- 
vide a microscopic foundation for the effective RMT, or 
alternatively, to derive Eqs. 0) and © from the qua- 
siclassical theory of Refs. yl |3- One might think that 
diffraction of a wave packet at the point contacts is the 
key ingredient which is presently missing from quasiclas- 
sics and which would eliminate the exponential te/td 
suppression factor from Eq. Q. However, our observa- 
tion of an exponential te /t^ suppression factor suggests 
otherwise: If diffraction at the edge of the point contacts 
were the dominant mechanism by which wave packets are 
split into partial waves, then the characteristic time scale 
for the suppression of quantum interference by dephasing 
would not be different from the mean dwell time td- 
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